Dynamics of Cosmological Perturbations in Multi-Speed Inflation 

Shi Pi^'H and Du Wang^'H 

^Department of Physics and State Key Laboratory of Nuclear Physics and Technology, 

Peking University, Beijing 100871, China 
(Dated: August 4, 2011) 



Abstract 

We study a multi-field inflationary theory with separable Lagrangian, which has different speed 

C^ I of sound for each field. We find that the fields always coupled at perturbative level through 

gravitational interaction. We show that if the coupling terms among the perturbation fields are 

weak enough, these fields can be treated as a combination of decoupled fields, which are similar to 



(N 






(N 



<^ . normal modes in coupled oscillation. By virtue of such fields, the curvature perturbation at the 

horizon-crossing can be calculated up to the leading order of slow variation parameters via SM 
formalism. Explicitly, we consider a model of multi-speed DBI inflation, and calculate the power 
spectrum in detail. The result depends on the ratio of different speeds of sound, and shows an 
apparent amplification when the ratio deviates from unity. 
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I. INTRODUCTION 

Inflationary cosmology has become the prevalent paradigm to understand the early stage 
of our universe, with its advantages of resolving the flatness, homogeneity and monopole 
problems [l|, |2|, and predicting a approximately scale- invariant primordial power spectrum 
consistent with current cosmological observations [3| very well. However, a single field 
inflation model often suffers from fine tuning problems on the parameters of its potential, 
such as the mass and the coupling constant. 

In recent years, people has noticed that, when a number of scalar fields are involved in 
the inflationary stage, they can relax many limits on the single scalar inflation model j^. 
Usually, these fields are able to work cooperatively to give an inflationary stage long enough, 
even none of them can sustain inflation separately. Models of this type have been considered 
later in Refs. |5|-l8|. The main results showed that both the e- folding number Me and the 
curvature perturbation TZ are approximately proportional to the number of the scalars A^. 
Later, the model of N-flation was proposed by Dimopoulos et al. [9|], which showed that a 
number of axions predicted by string theory can give rise to a radiatively stable inflation. 
This model has explored the possibility for an attractive embedding of multi-field inflation 
in string theory. 

Over the past several years, based on the recent developments in string theory, there 
have been many studies on its applications to the early universe in inflationary cosmology. 
An interesting inflation model, which has a non-canonical kinetic term inspired by string 
theory, was studied intensively in the literature. Due to a non-canonical kinetic term, the 
propagation of field fluctuations in this model is characterized by a sound speed parameter 
and the perturbations get freezed not on Hubble radius, but on the sound horizon instead. 
One specific realization of this type of models can be described by a Dirac-Born-Infeld-like 



(DBI) action [10|, lHj. Based on brane inflation [12|], the model with a single DBI field was 



investigated in detail|13Hl5|. which has explored a window of inflation models without flat 
potentials. In this model, a warping factor was applied to provide a speed limit which keeps 
the inflaton near the top of a potential even if the potential is steep. 

Motivated by the effective description of multiple D-brane dynamics in string theory 
16|, an interesting inflation model involving multiple sound speeds with each sound speed 



characterizing one field fluctuation was proposed [l7|, [l8|, |3J]. The authors of Refs.|l7l. llSl 
suggest this scenario can be realized by a number of general scalar fields with arbitrary 
kinetic forms, and these scalars have their own sound speeds respectively. Therefore, this 
model is dubbed as ^^Multi-Speed Inflation" (MSI). In this model, the propagations of field 
fluctuations are individual, and the usual conceptions in multi-field inflation models might 
be not suitable in this scenario. For example, in a usual generalized N-flation model, the 
length scale for perturbations being freezed takes the unique sound horizon; however, in 
our model it corresponds to the maximum sound horizon. It is worth emphasizing that the 
model of MSI is different from the usual DBI N-flation in which only multiple moduli fields 
are involved in one DBI action |l9l-|25|. but MSI is constructed by multiple Kessence type 



actions. For example, an explicit model of MSI is made of two DBI fields in Ref [17[, and its 
primordial perturbations including curvature and entropy modes and their non-Gaussianities 
were considered. 

In the current paper we extend the model into a much more generic case by relaxing the 
form of the Lagrangian. We find that even there are no coupling terms among the inflaton 
fields, their perturbation modes are coupled due to the nonlinear gravitational interaction. 



We also find that if tliese couplings are weak, the independent components of the fields can 
be treated as a combination of normal modes, by which the curvature perturbation at the 
horizon-crossing can be calculated in detail up to leading order of slow variation parameters 
via SM formalism. As a specific example, we consider a model of multi-speed inflation 
involving two DBI fields. We show that in the relativistic limit, the coupling between two 
fields is mainly reflected by the damping terms in the perturbation equations. The difference 
between the sound speeds of these two fields is able to generate a considerable amplification 
on the amplitude of the primordial curvature perturbation. This is greatly different from 
the usual analysis on the inflationary model constructed by two canonical fields(2j 
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The paper is organized as follows. In Section II, we study the cosmological perturbation 
theory of the MSI model in a generic case, then show that the infiaton fluctuations could 
be decoupled through a redefinition of fields at leading order under the assumption of weak 
couplings. In Section III, we focus on a model of MSI involving two fields, and study the 
detailed field transformation matrix to illustrate that the decoupling process of the field 
fluctuations is reliable under the weak coupling approximation. Specifically, we analyze a 
specific model constructed by two DBI fields and give its curvature perturbation. Section 
IV presents a summary and discussion. 

In the paper we take the normalization M^ = I/SttG = 1 and the sign of metric is 
adopted as (— , -|-, -|-, -|-) in the following. 



II. THE MSI MODEL 

An inflation model constructed with a single Kessence was originally proposed by [3 



and later its perturbation theory was developed in [3l|]. In the literature this type of model 



has been widely studied, and one of the most significant features is that there is an effective 
sound speed describing the propagation of the perturbations. In the model of MSI, for each 
k-essence field there is one sound speed correspondingly. Therefore, the field fluctuations in 
this model do not propagate synchronously due to different sound speed parameters, but get 
coupled because of gravitational interactions. In the current paper, our main interests focus 
on the effects of multiple sound speeds and mode-couplings at the level of linear perturba- 
tion. Before studying the perturbations, we first take an investigation on the background 
equations. 



A. The setup of background model 

The action is the summation of Hilbert-Einstein action and the action of A^ different 



fields minimally coupled with gravity [17 

S= I (fx 
with 
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which is the kinetic term of the infiaton field (pj. From now on the subscript of capital Latin 
letter I,J,... always denote the different infiaton fields. In this model, the Klein-Gordorn 



equations in a homogeneous background are 

-#^ = 0. (3) 




Here and thereafter, the subscript ", /" denotes the derivative with respect to the field 
0/, and ", J^/" denotes the derivative with respect to Xj, for simphcity. Assuming a fiat 
Friedman-Robertson- Walker background, the Friedman equations are given by 

3H' = p = J2i'^XjPi,x,-Pi), (4) 

I 

H = -\{p + Y.Pi) = -Y,XiPi,x,, (5) 

where H denotes the Hubble parameter a/ a at a given time. 

To get analytical results, it is helpful to define some parameters to characterize the kinetic 
behavior of the fields. In single-field canonical infiation, a set of slow-roll parameters are 
necessary. To ensure their smallness is not only the requirement of inflation to occur and 
endure for long enough, but also some simplification in analytical calculation. In a generic 



case, just as in [32], we define a similar set of slow- variation parameters. Recall that in the 
simplest single-field infiation with canonical kinetic energy and fiat potential, one can define 
two slow-roll parameters 



H 
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n = — ~ — r ~ . (7) 

' He H<p V^(0) ^ ^ 

Here the ~ symbol denotes the approximate equality which only holds in slow-variant envi- 
ronment and neglect the possible numerical factor. The second equality defines the slow-roll 
parameter via the kinetic side, whereas the third equality defines it in the potential view- 
point. If one only deal with the parameters under slow- roll condition, for instance in a 
model with canonical kinetic terms, these definitions are almost equal, and assuming one 
set is small can immediately assure the smallness of the other sets. This is just because the 
dynamic behavior simply relies on the shape of potential. But in a generic case, for instance 
in DBI or k-infiation model, the dependence may be much more complex and these different 
definitions are no more equivalent. Since they are all very important in parameterizing our 
final result, we have to define two sets of slow- variation parameters from different view- 
points. In general, we should replace the slow-variation parameters by some slow- variation 
matrixes that can reflect the correlations between different flelds. Since we are dealing with 
a separable Lagrangian, we only need different parameters for each fleld. This is equivalent 
to the case with a diagonal slow-variation matrix. As is mentioned above, one way to trace 
the motion is to deflne the parameter by each kinetic term of one speciflc fleld: 

^. = %^, (8) 



In DBI model the field may roll down the potential very fast, but the parameter e/ defined 
here is still small and slowly variant. The summation of such parameters gives the original 
version of them in canonical case: 

^ = -^ = E^- (10) 

V = Jf^ = ^'^i- (11) 

Note that the magnitude of e/ and e obeys e//e ~ 0{1/N) in general. To accomplish our 
calculation we need another version of slow-variation parameters defined in the potential 
aspect: 

e, - -i|. (12) 

h. ^ ^. (13) 

This is an analogue of e ~ iV' /VY, and for future convenience we take a squared root. It 
is not necessarily small in an infrared (IR) DBI model, in which the value of the infiaton 
field may be very small. So for convenience let us focus on the ultraviolet (UV) DBI model 
where (f) moves from UV side of the potential to the IR side. We should remind us that e] 
has the same order with ej in a model with canonical kinetic energy, but is very close to 1 
in multi-DBI model. We refer to it as a half-order slow variation parameter in the former 
case, and tried to preserve terms higher than ef. Besides, 

as another extension of ?7 ~ V" /V should be defined. According to the speed of sound in 
the single-field case, we can define the speed parameter of each field: 
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By this definition, the sound speed is the actually the effective propagation speed of the 
perturbations, which is different from the adiabatic sound speed which is defined as c; 
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P/p when the cosmic fluid is described by scalar fields (3J]. Notice that in a multi-Dirac-Born- 
Infeld model we have an extra relationship Pi,Xj = l/cs/ which will connect this parameter 
with e/. In general, Csi needs not to be the same for all the fields, which makes this model 
different from the ordinary multi-DBI model and may give us some new phenomenon that 
we are about to investigate. For an analysis of Cg as an inverse Lorentz vector in a single-field 
DBI inflation, see |29| . 

After having defined such parameters, the main quantities involving Pj can be expressed 



in terms of these parameters directly from the inverse relations above. For instance, 

Pi,x, = ^ej, (17) 

^— T(0(i-). P«) 

Pi,i = -SH^ei, (19) 

Pi,n = -3H\i. (20) 

Besides the P's, another important quantity is the time derivative of one field. By using the 
Klein-Gordon equation (^ we have 

^='-L(l-'-e'^- le^) . (21) 

2H ei\ 3 e/ 3 hj J ^ ' 

B. Generic perturbation analysis 

In this subsection we perform a generic analysis on the linear perturbations of the model 
introduced previously. Since we are working in the frame of a cosmological system, the 
metric perturbations ought to be included as well as the field fluctuations. We would like 
to expand the action to the second order by virtue of the Arnorwitt-Deser-Misner (ADM) 
formalism J35[. 

To start, the spacetime metric in the ADM formalism is written as, 

ds^ = -N'^df + hij {dx' + N'dt) {dx^ + N^dt) , (22) 

the action can reads 

S = - I dtd^xVhN I -i?(3) + J2^n^~ [ dtd^xVhN'^ [EijE"^ - E'^) (23) 

with Eij = (1/2) (/ijj — ViNj — VjNi). When expanding the action up to the second order, 
we can decompose the fields and Lagrangian multipliers as 

0,(t,x) = 0,(t)+50,(t,x), (24) 

A^ = 1 + a, (25) 

Ni = d,p. (26) 

First, the primary Hamiltonian and the momentum constrains gives two relations 

I 



Then we can make the perturbative expansion of the action up to second order by using the 
Lagrangian constraints: 
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dtd X— <j Pi,XiXi vj + Pi,xi S4>i 
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— Pl,Xj 0/50/ (Pj^Xj Vj + Pj,J 5(f) j) + Pi,ii I 
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vj = (j)i5(j)j - ( ^Pj^Xj (pi5(l)i 
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as the perturbation of Xi up to first order. This result is in consistency with the one derived 
in the MSI model 17| and see 36|, |37| for a general case. The Lagrangian fl28|) contains a 
lot of coupled terms as 6(f>i6(j)j and S(j)i6(f)j. These terms imply that the perturbations of 
inflaton fields depend on others during their evolution, and the equations of motion would 
be difficult to be solved directly. In this note we aim at making the transformation of the 
inflaton perturbations as to obtain approximately decoupled equations of motion for new 
flelds. In principle this is difficult, if not impossible to realize. We will see in which case 
can we do such transformations and see what will happen by such a method. This will 
require some constraints on the slow- variation parameters. We will see this in detail for a 
double-field toy model. 

After utilizing the equations of motion for background and taking some integration by 
parts, one can write the Lagrangian in a more compact form: 
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dtd xa 
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[Pi,Xi +Pi,XiXi ^j)S<l)j ^ 



ij<- 



Pi,xj di6(j)idi6(j)i 



(30) 



with cAjj and ^u are time-dependent damping and mass terms, 
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(32) 



This result is consistent with that obtained in multi- field infiation 23|, |38|, |39| . To distinguish 
the diagonal and off-diagonal parts of the matrices above, we define 



^ij = U6ij + 7V/J, ^1.1 = M6ij + Mij. 



(33) 



After assuming that all the slow- varying parameters defined in Sec. |T] are small, we have 
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Here, note that all these quantities are of order 0{e) if c^j ~ C^(l)- I^^ such a case, as we 
will see below, the damping term has no correlations. For ordinary terms we have only 
preserved the leading order to slow- variation parameters, but if a term involves c~j^, we 
preserve one higher order. This is to ensure the consistency when c^j ~ 0{1) fails, as in 
a model with multiple DBI actions. Under the condition when c^j ~ 0(e), all the terms 
should be preserved. 

To investigate further we can derive the equations of motion for the fields 
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Note there is no summation over index "/" even if it repeats in one term. Here and through- 
out the paper Einstein's summation rule is not applied for the index of fields (but still valid 
for repetition of spatial sub /superscripts) and all the summations over capital Latin letters 
"/, J, ..." are written explicitly. In a generic case, we can rewrite the equations of motion as 
follows, 



3H{1 + Ki)5<Pi + [^ + m\ 64>i + HY,iiMj + $^m,j50j = 0. 



(39) 
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mi 



mij 



Ki 



:iJ 



"s/ 



Pi 



■M, 



,Xi 



"si 



Pi, 



Xi 



Y.^ij^ 



2e 



-l^Sl + -7]i - -7]iC,j - — 1 



VI 



e/ 



2ei—c^j ( -^ 



ei 



-sj 



hi 

ei 



"slJy 



(40) 
(41) 

(42) 
(43) 



which are also small in inflation. We see there are two different couplings in the damping 
and effective mass term respectively. Besides, they both depend on time via the parameters. 
In principle, this kind of equation is difficult to solve, if not impossible. Of course we see that 
the coupling terms between different fields are all very weak if we admire the slow- variation 
parameters to be small, which may maintein the final solutions almost the same as those of 
N independently evolving fields. Further more, two special cases which have specific limits 
of sounds of speed are most interesting. One is that all fields (pj have the canonical kinetic 
energy, with Cg/ — )■ 1. This is a decoupled multi-slow- roll model, which is called assisted 
inflation by [4] and received much focus[5|-l8|. In such a case, ^u ~ 0{e^) and therefore can 
be totally neglected in the action, while mu ~ 0{e). We notice that only the correlations 
in mass term are preserved, which are easy to manipulate. On the other hand, when all 
the fields are DBI-type, as is discussed in [l8|, Cg/ is very small in the relativistic limit. For 
simplicity, we can suppose that the sound speeds are as small as e, or even smaller. We will 
shortly confirm that, the latter case is naturally suitable for analytical calculation, while for 
former case one should impose further assumptions on the difference of Cg's. But whatever 
case it is, we see ^u ~ 0{e) which is the main part of coupling, and mu < 0{e'/'^) thus can 
be neglected. 

One possible method is to consider the smallness of couplings, and rotate the perturbative 
fields in field-space into a instantaneous orthogonal basis which can generate the adiabatic 
perturbation which is in the direction of the trajectory of fields in background, along with 
N — 1 entropy perturbations perpendicular to the direction of field trajectory. After doing 
that, if we suppose the coupling terms which are proportional to the derivatives of the ro- 
tating angles with respect to time of the trajectory in field-space is small at the moment of 
horizon crossing, we can neglect the coupling and treat the adiabatic and entropy perturba- 
tions independently as free streaming quantum fields. After some proper quantization, all 
the discussions are similar to a A^ decoupled multi-field theory [231. l38l . |39|. Next if we want 



to consider the effect of couplings that have been omitted, one can treat such couplings 
as interacting vertices to the quantum fields, and see how can the fields transfer to each 
other in the quantum level [40|]. Actually, if the fields have different speeds of sound, it is 
difficult to do such rotations as to make the adiabatic perturbation laid on the direction of 
trajectory. This is because the angle we need to take our rotation is different from that of 
the slope of trajectory, since sound speeds will appear in the transformation and make the 
"rotation" anisotropic. Further more, different sound speeds correspond to an anisotropic 
rescaling in the rotation of basis. We will come back and explain this assert in a two-field 
case by illustration later. 

If we abandon the attempts to project the field into adiabatic perturbation and entropy 
perturbations, we may find another route to get our results. A reasonable method is to 
transform the fields in order to eliminate the couplings, if possible, and then treat the 
decoupled perturbations as free quantum fields, with a final result of curvature perturbation 
TZ gained by dM formalism. Before doing this, let us digress to take a look at the solution if 
we discard all the couplings and see the free fields limit of our theory. The result obtained 
here is also useful for future convenience. Under this situation, the equations contains no 
coupling, while the only remaining task is to find the solution of each independent field, and 
do quantization in ultraviolet limit. We can consider the equations of motion as 

5(t)j + 3/7(1 + Ki)5(l), + ( — ^ + m] \ 5(t)i = 0, (44) 

where the parameter k and m have been defined in P2l) and (j^Ol) . This is a little different 



from the standard form of single-field inflation because of an additional term in the damping 
effect. We consider the fields moving in an inflationary quasi-de Sitter phase, which means 



the comoving time has a simple form [41 



Hal 
And we can postulating a rescaling on the fields, 
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Note that, our rescaling contains an extra power of /3 which is a small quan tity of order 0{e) 
comparing with the standard definition of Mukhanov- Sasaki variable 42| due to the small 
shift in the diagonal damping term of the Lagrangian. This extra power can be determined 
by requiring the coefficients of x'l term in the equation of motion vanish. The Mukhanov- 
Sasaki variable, as in the standard process before taking the quantization of the fields, is 
actually defined via the normalizer before the kinetic terms in the action f l5U]) . Therefore, 
we can convert the equation into a standard form which is familiar to that in single field 
theory: 



x'; + 



k Cgj 



1 



9 3 



2IP 



XI = 0. 



(49) 



This is just the Bessel equation of order uj which will be determined below. The solution to 
this equation is a linear combination of two independent Bessel functions. However, to have 
the appropriate approximation behavior which approaches the planar wave in Bunch-Davies 
vacuum with correct normalization at early time r — )■ — oo J43|, we need to choose the proper 
coefficients of the Bessel functions, which corresponds to such a state in quantum theory 



that it can minimize the energy density. One possible result is to choose J4l|, |42 



Xi 
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where Hu/ is the Hankel function of the first type, with 
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(51) 
(52) 



as its order. Note that, o"/ is a small quantity of order 0{e), which is irrelevant of the value 
of Csi in a multi-DBI model when Cgj ~ 0{e). Therefore, the amplitude of perturbation to 
the original field 0/ is 



-,1^1-2 _ 



r(^'/) ( c^i^fWi 



r(3/2) 



aei 



2(l+ft) 



-kc^jT) 
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(53) 
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Again we use the fact that both vj and k/ are very small, thus we get the power spectrum 
of each field 



V, 



i^ — [1 + (/5/ - 2^/) 1^ 2 - a,^ (3/2)] ( 5 



\2'K J Csl 
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e/ 



l-2e + /3/ln2 + ffj(7 
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2) + /37 
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N2(l-a,) 



(54) 



where the second approximate equality holds at the moment when the wavelength of the 
mode considered exits the sound horizon, i.e. kc^i = Ha. ip is the digamma function and 
relates to Euler-Mascheroni constant 7 as ip{3/2) = —7 + 2 — 2 In 2. The dependence on the 
speed of sound is superficially different from that of single-field DBI inflation. But when we 
take an extra relationship H = — Y2jXi/csi which holds in DBI-type action, a new relation 
between the parameters Cg/, e/ and e/, holds as to leading order. 



Csl- 



(55) 



Therefore the dependence on Cgi in ( l54l) is canceled, and (IMI) will be coincident up to leading 



order with the result in (32|. By virtue of SJV formalism 4^, |45 . 

/ ij 

where A/'e is the local e-folding number along a trajectory from a spatially fiat slice at a 
moment t* soon after the relevant scale has passed outside the horizon during inflation, to 
a uniform- density slice at another moment tc after complete reheating when ( has become 
a constant. So the power spectrum curvature of perturbation ( can be calculated by |44| . 
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H^e,] 



where iVe,/ can be determined via 
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I ^^ I 

in a general case. 

Here for simplicity, we should consider the case when the A^ fields moves almost in 
the same manner. This is proved by |j] to be a late-time attractor for some appropriate 
initial conditions. Later on we will demonstrate our analytical calculation for a system with 
different sound speeds is valid when e/ — ej ~ C'(e^), which implies a field configuration like 
this. Under such assumptions, to order 0{e) one have 



dAfe 



1 ej 



T 1 ^/ 2 e/ 

1 H — e e-— 

3 e/ 3 hj 



d(t)i 2N ei 

in a generic case. Therefore, together with f l55|) we get 

hV 1 ^ 1 /i/^'^ 
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The superscript denotes this power spectrum is the one for a theory neglecting couphngs. 
We can see if all the fields have the same speed of sound, and e/ ~ e/iV, ( 160|) gives the 



results in multifield DBI inflation [23[ which contains A^ identifying speeds of sound up to 
first order of slow variation parameters. 
From (1541) . the scalar spectral index ris is 

ns7 - 1 = -2/3, = -2e(l-^- ^) (1 - c^,) + 2r/,(l - c^,) - 2hi. (61) 

And the contribution to spectral index of ( mainly comes from the largest one of rig/'s, which 
is just the power spectuum of curvature perturbation. 



C. Define New Fields 

To consider the coupling terms as a perturbations to equations of motion fl39|) . we could 
take a redefinition of infiaton fields through a linear transformation which, up to leading 
order in slow- variation parameters, can give a set of N decoupled equations of motion, with 
different coefficients appearing as new speeds of sound and new effective mass terms. In 
general, this is difficult to be realized. Fortunately, we are dealing with a theory that is just 
perturbed from a decoupled one. As the A^ fields have tendency to move together, they have 
an even smaller difference in the parameters. Intuitively the change in eigenvalue should 
be small if the damping and mass matrix is not far apart from a diagonal one. Actually 
Hoffman- Wielandt theorem in matrix analysis tells us that the smallness of "perturbation" 
to a diagonal matrix will cause small variation to the eigenvalues. To make it clear let 
us denote A/ as the eigenvalues of a diagonal matrix corresponding to a theory with no 
couplings, whose diagonal entries are just the coefficients before each 6(f) term in ( 139|) . and 
A/ to be eigenvalues of the "perturbed" matrices corresponds to the coupled coefficients 



Then, there exists a permutation a (J) for the new eigenvalues A that satisfies (46 

1/2 



■ TV 

E 

7=1 



A. 



<\\6M \\2r-0{Ne). (62) 



Here || 6M ||2 is the Euclidean norm of a perturbation matrix 6M which under our circum- 
stance has the magnitude of order 0{e). Thus, roughly speaking, the variation in eigenvalue 
is a small quantity of order less than 0{e), if A^ is not very large. This fact ensures our 
validity of our upcoming discussions on estimating the eigenvalues after field redefinition. 

Let us analyze the coupled equations under the fact that all the coupling terms are small 
quantities of order C(e). In this case we expect the solutions will differ from the decoupled 
ones very little. Later on we will see that, even under such case, different sound speeds will 
bring a fruitful physical imprint. If we impose a transformation in field space such that 



6h = Y,RTfi, (63) 



where R is a non-singular transformation matrix. If we require the transformation to be a 
representation of SO{N) we can define a new basis of perturbations which can be interpreted 
as adiabatic/entropy directions, this matrix should be set orthogonal. But here it is not 
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possible in our case. Our purpose of transforming the fields is to decouple. And unless for 
some special cases, the new fields has nothing to do with adiabatic/entropy perturbations. 
Actually, if all the fields have the same speed of sound, the fields redefinition (163|) could 
be an orthogonal transformation of SO{N) group. And one of the new basis would lay on 
the direction of tangent line of trajectory in field space, which just duplicate the result in 
28| . However, for fields with different speeds of sound, we do not have such clear geometric 



explanation. In such a case, we will see that the transformation is even only orthogonal up 
to leading order 0{1). Later on we will see how we can explain the new transformation in 
a field space. 

To make our discussion clear we could require R is also slow- variant, which means M ~ 
0{e). This is reasonable in a quasi-de Sitter background, and its validity will be confirmed 
after we get the concrete form of M. Thus, one can write the time derivative of new field 
transformed from M.S(b as 



bj ^ R^^Pm + IRT^m. 



(64) 
(65) 



Substitute into the equations of motion fl38p . we have the equations for new variable </?. Left 
multiply an inverse transformation (i?"^)^ and take the summation over index J, we get 



m,I,J 
m,I,J 

+ mjjRJ + H^ijR 



3Hil + Kj) + 2RY)5u + H^jjR 



^r 



^i?7 + mJRy + 3/7(1 + Ki)R'y ) Su 



^r, 



0. 



(66) 



If we need to decouple these equations into the form of ( l44l) and get a set of free propa- 
gating waves in ultraviolet, we should find a way to diagonalize the two coefficients in the 
square brackets before and ip for each equation at the same time. This is the simultaneous 
diagonalization of two matrices, and is in principle impossible unless the two matrices are 
commutative. The coefficient matrices does not commute each other since E,ij is even asym- 
metric. All the diagonal matrices are commutative. So the only remaining matrices that 
may be non- commutative are those with m/j and ^/j as their elements. Since their main 
part are diagonal elements that surely commutes, we will justify that the complete matrices 
with a little variation in off-diagonal parts are believed to simultaneously diagonalizable up 
to 0{e). If we can find a way to do this, then we can write 



E 



2{R-y„RT + 3H{R-y„K,RT + H^(R-XCuR 



E 



3H(R-X(1 + ^,)RT + -(R-'y„ci,Rf 



+(i?"i)>?i?r + j2iR~'yn^hR 



3HKJl^, 



(67) 



-Cn + M„ C- (68 
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Here we have neglect a ^M term since M is of order 0{e). Under such transformation the 
equations of motion are completely decoupled, 

^n + SH{1 + Kn) + [KCsn + MA ^n = 0. (69) 

Now we use the Hoffman- Wielandt theorem which ensures us that the solution to such 
eigenvalue equation of a perturbed matrix is a small variation from the eigenvalue before. 
According to the fundamental theorem of algebra, there are N sets of such solutions to 
dnZ]) and (jHED, each has a new speed of sound C^. The Cg^'s are the speeds of propagation 
of the new fields, akin to those of normal modes in coupled oscillation system or phonon 
representation in statistical system. To investigate whether the new representation will 
change qualitatively the physics of an decoupled system is an interesting question. For 
instance when Cs/'s are all close to 1, can we find a mode that has a very small speed of 
sound Csn that can possibly amplify the power spectrum and non-Gaussianity? We should 
emphasize that although this is possible, it will not appear in an weakly coupled system 
we are just dealing with since the reason we've mentioned above. On the contrary, if the 
couplings are indeed strong, the physics could be much more different. For example, in 
hybrid inflation 47| where the two fields have complicated coupling which is not small at the 
waterfall stage of inflation, there are plentiful physical phenomena far beyond our results 
here |48l451|. Even in a weakly coupled system we are considering, analytical calculation is 



difficult to do. But after we can diagonalize the damping and mass matrices and take a set 
of new fields </?„ which can simplify the equations of motion fl38|) to f l69|) . the remaining task 
is very simple, just parallel to the discussion under dH]) since the form of the equation is 
the same, only with a different definition of parameters Kn-, Cgn and Mn. The most difficult 
part of our method is to find the transformation M. In the next section we will calculate in 
detail for an example of two-field infiation and get a taste of such scheme. 



III. TWO-FIELD CASE 

To calculate the curvature perturbation at horizon-crossing in detail we need the quan- 
tization of the fields in ultraviolet, which requires a set of independent wave equations. As 
we mentioned above, to simultaneously diagonalize so many matrices in f E7|l and fl^Sjl seems 
to be an impossible task. But after we notice all the off-diagonal entries are very small and 
superadd an assumption that the diagonal entries are almost the same up to 0{e), we can 
decouple the equations in a natural manner. In a two-field model, the fields which drive 
infiation are 0i and 02, and a new definition of fields provides two equations for the eigen- 
vector of the coefficient matrix, ( 167|) (168|) . We have seen that for a generic value of Cg/ the 
off-diagonal damping and mass term will coexist. But for a model with canonical kinetic 



energy, Cgi — >■ 1 and there is only coupling in mass term|52H54|. On the other hand if Cg/ is 
small as in the multi-DBI model, the coupling in damping term survives. Let us work in the 
multi-DBI case, leaving aside the matrix nijj which has been proved to be up to (9(e^/^), as 
has been discussed in the paragraph under fHT]) . Therefore the only remaining matrix to be 
diagonalized is ^u, which has the eigenvalues 



■^21 'C22 — ^ 
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0. (70) 



With the eigenvalues 

AW = 2(ei + e2), A^^^ = 0. 

These eigenvalues corresponds to a transformation matrix 



£2 
ei 



-1/2 



/ 



V 



€261 C 2 

eie24 



e2C,i \ 
eic,^2 



/ 



(71) 



(72) 



Actually the eigenvector can only determine this matrix up to two normalization constants. 
Here we adopt a normalization such that it will preserve the determinant to be unity. This is 
not important in the quantization process, but will affect the calculation of power spectrum 
later. We know the similarity transformation M will diagonalize the matrix with elements ^/j. 
Next we will show that it can also diagonalize the diagonal matrices with different diagonal 
entries, say K = diag(fi;i, /t;2), under some reasonable assumptions about the model. When 
the similarity matrix R operates on K, 



(73) 



1 




eiKi + 62^2 

■7 ¥(K9 — Ki 


) ei Ko + eo 


-Ki)^ 


ei + 62 



So a diagonal matrix which has different diagonal entries can remain diagonal after the 
similarity transformation Mif 

/ti - K2 = O(e') (74) 

with K/'s and e/'s in the matrix are still of order e. We will see, that for a specific multi-DBI 
model, the assumption that ni — K2 as well as Si — S2 and r/i — 772 are of order O^e"^) is 
reasonable, since the difference between two small positive quantities of order 0{e) must 
be even smaller, if the fields have the tendency to move together under an appropriate 
initial condition. Besides, we will suppose that the differences between other parameters 
like ml — m\ and (^-^ — 0^2, are also as small as 0{e^) for future convenience. Note, that this 
do not implies ei — £2 ~ 0{^^) which is a much more strict constraint on the evolution of 
fields. We will come back to this issue and make some necessary justification to a specific 
two-field-DBI model later. At this moment, we just emphasize that this assumption preserve 
the diagonal property of the diagonal matrices in fl67|) and fl68l) involving k/, c^j and ?n| 
automatically up to 0{e^) under a similarity transformation by any matrix M. 

Next let us check whether this matrix can make M^^M diagonal. The entries of the matrix 



are 



He2 



[)-1t 



)-1t 



ill 



;i2 



;2i 



i22 



ei + £2 



Hei 620. 



si 



ei + €2 eic; 



1 

2" 

^2)] 



•s2 



He2 eic^2 
ei + €2 eac^^ 



{hi - /la) - -(?7i - r/2) - 2(si - S2) 
[(/ii-/i2)-2(si 
[{hi - hi) - {rji - 



V2) 



2{Si-S2)]. 



(75) 
(76) 
(77) 



We see all the entries are proportional to the differences among some slow-variation param- 
eters of the two fields, which are of order 0{e^) as we mentioned before. Then we can assert 
that up to 0{e) the time derivative of M multiplied by its inverse in ( 167|) is negligible. 
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54>2 


V^i ^ 


' 


1 ^ 






/ , 5(6] 




/ 




01 




(a) 



(b) 



FIG. 1: The transformation M is depicted by a transformation in field space of the perturbations, 
(a) denotes the trivial case when ei = e2, ei = 62 and Cgi = 0^2- We see for a totally symmetric 
configuration M poses an S0(2) rotation by an angle '7r/4, and 0i is just the adiabatic perturbation 
in usual case, (b) shows a more interesting result when ei = (l/-v/2)e2, ei = (l/\/2)e2 and 
Csi = \/2cs2. Note this is to ensure 4>\/(l)^ = eiCs2/(e2Csi) = 1 such that field are moving along 
a diagonal trajectory in fields space. We see that M represents a transformation that rotates the 
field and rescales it at the same time with different rotating angle and anisotropic rescaling factor 
for two fields, i.e. the red arrow denoting the basis of new fields are not orthogonal, and the new 
fieds v? can not preserve the unitarity. 



We will pause for a while and see the possible geometric explanation of transformation M. 
We have seen that it is only orthogonal up to C(e) under the assumption that ei —€2 ~ C'(e^)- 
When Csi = 0^2 as is the case of ordinary assisted inflation, one have that the transformation 
M is just rotating the basis by 7r/4 and ip2 is just the adiabatic perturbation with ^pi the 
entropy perturbation which is negligible, as is depicted in FigJII^a). A more interesting case 
is the theory with different speeds of sound, as in (b). 

After the transformation, one can find that two equations (!38|) are decoupled, which is 

just iSD, 



(pn + 3i7 (1 + i^„) + ( ^Csn + Mr 



^r, 







with 



Ki 



eiKi + e2K2 
ei + £2 

p2 



ei + 62 



^iCsl 



e2Cs2 



Ml 



ei + e2 
eim\ + e2m2 



K, 



'-^s2 



Mo 



€2^1 + 61^2 



ei + 62 






^ ^2 



ei + £2 
e2mi + eim2 

61 + 62 



(78) 
(79) 
(80) 



ei + 62 

To go further let us study the change of the normalizer zj under the similarity transformation 
by M. Note that, for example, when a normalizer matrix Z, which is also diagonal, is 
multiplied by another matrix, say X, the product is transformed under M by 

R-iZXM = M~^ZMM~^XM, (81) 
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which is just multiplying each diagonal entries after the similarity transformation respec- 
tively. Denote by Zi the diagonal entries of ]R~^ZM, 






aij) Cs2eiy^+ Csie2A/e^ 



a(T) Csie2e2ei 



CslCs2(ei + £2 

1/2 



Cs2eieie2 



1/2 



V^ 



CslCs2(ei + £2 J 



After doing this, we can turn to the same process after equation 
field variables 



B2 



3 _ Cs2eiv/ei Ai +Csie2A/e2A2 

2 Cs2ei^+ Csle2^/e2 

o -1/2 A , -1/2 A 

3 _ Cs2eieie2 Ai +Csie2e2ei A2 



-i^. 



Csie2e2e 



'1/2 



Cs2eieie2 



1/2 



-^{r]-rin) + hn 



(82) 

(83) 
Then we define new 

(84) 
(85) 

(86) 
(87) 



which satisfy the equations of motion 



Xn ' 



.2^2 

''STI 



rC L^a 



9 3 1 

1 + -K„ + -e + ri A, 

4 2' 2 



M^ 



Ml 



2H^ H^ 



v) 



XI = 0. 



Now the equations have been decoupled. For each equation, the solution is a Hankel function 
with order 3/2 — ^„, where 



__ _2 _8 2Ml 
Finally, we obtain the power spectra of the new fields ipi as follows. 



(89) 



V 



Hy 2 



ka. 



2?i 



^^2(^1 + ^2)^ 



27ry Csi \ Ha J {Cs2eiy/e^ + Csie2y/e^) 



1 - 2e - 2(^1 (ln2 + V' (3/2)) + (3A'i - 2Ai) In 



gsigs2(ei + e2) 
Cs2ei v^ + Csle2^/e^_ 



(90) 



n 



^2 



H 



ka 



s2 



2?2 



Cs'i(ei + e2)' 



27r; ^2 V ^a 7 (^162626-'/' + ^261616-'/')^ 



1 - 2e - 2^2 (In 2 + ^ (3/2)) + {3K2 - 2 A2) In ■ 



CsiCs2(ei + £2) 



^81626261 ^'^^ + Cs2eieie2 



1/2 



(91) 



The power spectra above are calculated at the moment when the wavelength exceeds the 
sound horizon, respectively. When one field, say cpi have exceeded the sound horizon, the 
power spectrum of ipi is frozen, and so the curvature perturbation calculated later should 
involve the value oiV^^ calculated at t*i which satisfies Ha{t^i) = kCsi- The power spectrum 
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of curvature perturbation should be calculated at the moment when the wavelength has 
exceeded all the sound horizons corresponding to different (f^s. The power spectra are 
nearly scale- invariant, with — 2^„ as their spectral index. The spectral index of curvature 
perturbation ( is thus — 2<;"„ which has the largest absolute value. 

Using 6^/ formalism, we can calculate the curvature perturbation after the wavelength 
being stretched out of the largest sound horizon. Note that to leading order 

C = 6K = Y,^e,lS<Pl = Y,^e,lRTV>m. (92) 

/ I 

There is no correlation between different y^'s since they are independent kinetic modes as 
we defined. After using the additional relation (1551) we can substitute e„ in the expression 
by A/2e^/csn- Therefore we can simplify the final result 

n = (^e,l +^A4 J ) V^, + f ^A4,? +A4 J ) V^,. (93) 

The derivatives of the e-folding number can be calculated in a specific model in detail, for 
instance in a model with standard kinetic energy and separable potential|55|. Since we only 
want to emphasize the effects caused by different sound speeds, again we use the assumption 
that both fields are moving in the same manner, i.e. in the diagonal trajectory in field space 
as in (b) of FigJH and f l59p holds. This kind of motion is nearly adiabatic, which will always 
suppress the isocurvature perturbation inside the sound horizon. The final result is a little 
lengthy and instead of writing it here, we would rather focus on the dependence on different 
sound speeds by calculating Vc^/V)- . By defining Csi = xcs2 one gets 



x^ 



CL-^. CL-YT^- ^^^^ 



Now the power spectrum degenerates to 




(95) 

which now just contains quantities up to leading order, and is in coincidence with fj60|) when 
a; = 1, i.e. Cgi = Cs2- This normalized power spectrum seems to be divergent when a; — t- or 
X — )■ oo. However, since we require c^^ — c^j ~ 0{e'^), the corresponding physical region is 
the neighborhood of a; ~ 0{1). We depict the dependence on x in Figj2l 

We can see from Figj2]that, there is a minimum of P^ when Csi = 0.85cs2, where P^ = 

0.90P} . But the most interesting aspect is that if one of the speed of sound is a few times 
larger than the other (with themselves still as small as 0{e^) so as to satisfy our condition 
required before), the power spectrum of curvature perturbation increases magnificently. We 
see from Fig 12] and (!95l) that if a; > 1, the contribution from ip2 exceeds that from ipi, and 
vis visa if x < 1. Therefore we have encountered an interesting feature different from our 
intuition out of a theory without coupling, that the field with smaller speed of sound will 

18 




Csl 
Cs2 



FIG. 2: The dependence on Cg of the power spectrum. Both the power spectra of if are normahzed 
by {H/2'k)'^ . and Vc^ is normahzed by the un-coupled one (I60p . There is severe divergence of Vq 
when one of the sound speed vanish, but the only vahd domain under our assumption is the case 
when Csl ~ Cs2, i.e. in the neighborhood of x ~ 1. We see the amphfication by the reciprocal sound 
speed is much larger than that in a trivial x = 1 case of ordinary unique-speed model. 



not in general contribute more to the curvature perturbation. For instance, if c^i > Cs2, it 
is the power spectrum from auxiliary field (p2 which contributes more to the final curvature 
perturbation, not 6(f)i. We also see that if Csi/cs2 deviates from 1, the power spectrum of 
curvature perturbation is amplified significantly, compared with the one without couplings. 
This just tells us that the effect given by the coupling terms may be large and can not be 
neglect in previous calculations. Since the power spectrum of curvature perturbation can 
be observed to rather high accuracy, this result actually tells us that we can loosen the 
constraints on the slow- variation parameters of individual infiaton even further. 

IV. CONCLUSIONS AND DISCUSSION 

In the frame of inflationary cosmology involving multiple infiaton fields, it is possible to 
allow these fields possess their own sound speed parameters which characterize the prop- 
agation of field fluctuations correspondingly. It is important to study how these modes 
are related to the cosmological perturbation we may observe in CMB experiments. In the 
present paper we analyzed the dynamics of field fluctuations in a general scenario of MSI. 
Assuming no direct coupling terms among the infiaton fields, we noticed that these pertur- 
bation modes are generically coupled both through their effective mass terms and cosmic 
damping terms. Fortunately, we found that if these couplings are weak, the perturbation 
theory can be treated as a combination of normal modes, by which the curvature perturba- 
tion at the horizon-crossing can be calculated in detail up to leading order of slow- variation 
parameters via SN' formalism for a specific model. 

Specifically we studied a model of multi-speed inflation involving two DBI fields with their 
sound speeds being small and slow- varying. Our analysis showed that in the relativistic limit, 
the coupling between two field fluctuations is mainly contributed by the damping terms in 
their perturbation equations. We further studied the curvature perturbation in this model. 
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and verified their primordial power spectra are nearly scale- invariant. This has interesting 



implications to the curvaton scenario |56|. which suggests a light field in infiationary phase 



could seed entropy perturbation and convert it into curvature perturbation at the end of 



infiation[56H60| . Namely, as shown in |62H64| . a light DBI field can realize the curvaton 



scenario and generate sizable non-Gaussianities of local and equilateral types after suitable 
fine-tuning on the model parameters. Our study also showed that, when the sound speed 
parameters for the infiaton fields are not equal, the perturbation modes would get frozen at 
different sound horizons. An important phenomenon is that the curvature perturbation of 
MSI could obtain an enhancement which depends on the ratio of the sound speeds. 

The scenario of MSI is a very important branch of various models of multiple field in- 

fiation. In recent years, the infiationary models involving multiple field components have 

been studied extensively in the literature, namely, the analysis on the dynamics of N-fiation 

65| : multiple field infiation with particle decays J66|; the N-flation model in the frame of 



braneworld |67[; the model of staggered Nfiation in the stringy landscape [68|; djTiamics of 
the model of coupled N-fiation with canonical fields l69j; the infiation model constructed 



by a DBI field coupled to radiation (70[; and see |71l . |72[ | for comprehensive reviews on this 
field. It has been well realized that, for the model of multi-field infiation, when the sound 
speeds for the fields are the same, the field fiuctuations are able to be rotated orthogonally 
to a basis in which the curvature and entropy fiuctuations decouple explicitly, at least up to 
leading order. However, when we relax the scenario by allowing the sound speed parameters 
are different, our analysis implied that the rotation of the field fiuctuations are no longer 
orthogonal but we still are able to find a transformation to decouple the fields under some 
conditions and obtain the curvature perturbation via dM formalism. 
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